Correlation effects in disordered conductors with spin accumulation 
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We consider the effect of electron-electron interaction on the density of states of disordered param- 
agnetic conductor in the presence of spin accumulation and magnetic field. We show that interaction 
correction to electron density of states of the paramagnet may exhibit singularities at energies cor- 
responding to the difference between chemical potentials of electrons with opposite spins. We also 
discuss correlation effects on conductivity in metallic as well as in hopping regimes and show that 
spin accumulation leads to the negative magnetoconductivity. 

PACS numbers: 72.25. Hg, 71.20.-b, 73.23.-b 



INTRODUCTION 

The realization of large spin accumulation, that can 
be described by introducing the quasi-chemical poten- 
tials fj,±5[j,/2 for spin- up and spin-down electrons in the 
system possessing long spin relaxation time, is impor- 
tant for many proposals in spintronics. In particular, the 
non-equilibrium spin polarization might be created by il- 
luminating the sample with circularly polarized light [l[ , 
or can be achieved in conductors placed in contact to the 
ferromagnet via the spin injection mechanism 0, Q , see 
for a review [1, The spin injection and detection were 
investigated in many works and for different types of ma- 
terials such as superconductors organic polymers 
0, graphene [13]. Recently there has been a progress in 
achieving of large spin accumulation by means of electri- 
cal spin injection in semiconductors such as GaAs, Si, Ge 
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Motivated by experiments on spin accumulation, we 
focus on the non-equilibrium correlation effects in con- 
ductivity and density of states (DOS) in the presence of 
spin accumulation. 

Quantum corrections to the transport and thermody- 
namical properties of disordered metallic conductors have 
been a subject of both theoretical and experimental stud- 
ies. Electron-electron interaction in disordered conduc- 
tors results in the singularities of electron DOS at the 
Fermi level and positive magnetoresistivity, for a review 
see 
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For example, the energy dependence of the interac- 
tion correction to the electron DOS at zero temper- 
ature has a logarithmic singularity in two dimensions 
5u(e) oc In |er|, where r is the electron mean free time. In 
addition, external magnetic field B due to Zeeman split- 
ting leads to singularities shifted from the Fermi energy 
by the amount ±fl z as 5v{e) oc ln|(e 2 — 17 2 )t 2 |, where 
il z = \g/j,BB\sign(B) and [Ib is the Bohr magneton. We 
define Q z in such a way that it can be of either sign, 
depending on the direction of magnetic field. The inter- 
action correction to the conductivity decreases in mag- 
netic field in the limit of weak spin relaxation as [15| 



8a ~ — Aj/ln|ri z r|, where A is the electron-electron inter- 
action constant and v is DOS at the Fermi level per one 
spin. More detailed discussion can be found in 3, H|. 

Correlation effects related to the spin degree of free- 
dom can also play a significant role in magnetoresistance 
in the hopping conductivity regime. If the Coulomb re- 
pulsion between electrons permits the onsite double oc- 
cupancy, then the probability of certain transitions under 
the Zeeman splitting decreases leading to positive mag- 
netoresistivity lti- 



In the following sections we consider correlation effects 
in DOS and conductivity in case of finite spin accumu- 
lation (5/j,. We essentially rely on the calculations done 
for the equilibrium case, therefore we discuss only the 
features that distinguish the non-equilibrium case. 



DENSITY OF STATES IN METALLIC REGION 

We consider disordered metallic conductor in the pres- 
ence of spin accumulation and magnetic field. The spin 
accumulation can be obtained either by spin injection or 
optical orientation methods. We assume that the spin 
relaxation time in the system is much longer than the 
energy relaxation time such that we can treat the sys- 
tem in the energy equilibrium regime while having the 
non-equilibrium spin polarization. Within this assump- 
tion the interaction correction to the one particle DOS 
in the non-equilibrium state generated by spin accumu- 
lation can be treated in the similar way done for the 
equilibrium case 



14j . We need to take into account the 



non-equilibrium exchange splitting in the definition of re- 
tarded and advanced Green functions and to modify the 
Fermi distribution function of electrons by introducing 
the chemical potential shifts of electrons with spin-up 
and spin-down. 

Interaction correction to the one particle DOS is de- 
termined by the advanced component of Green function 



8v a {e) = -Im [ -^-SG A (p,£,a 

7T J (277)^ 



(1) 



where a = ± defines the spin direction and we accept 
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where Vb(q) is the Coulomb potential and IIo(w,q) 



Dq 2 /(Dq 2 



iuj 



Function F(oj, a) is determined by the 



Keldysh component of the non-equilibrium Green func- 
tion, which after averaging over random potential be- 
comes 



G K (P, ui, a) = —G R (p, uj, a) F (ui, a) G A (p, ui, a) (4) 

IT 



FIG. 1. Diagram for the calculation of the density of states. 



the units Ti = 1 . Let us consider the exchange part of the 
interaction correction shown in Fig. ([1}, where diffusion 
contribution to the self-energy part is given as 

S A (p, e, a) = G R (p, 0, a) — ^ J ^ x 

x [F(e, a) - F(u, a)} D% (u - e, q) V A (cj - e, q) (2) 

here D (w,q) = (-Dq 2 — iuj)~ x is the diffusion propaga- 
tor, D is the diffusion coefficient, V is the volume of the 
system and screened Coulomb potential is given by the 
following expression 



V A (uJ,q) = 



Vo(q) 



l + 2i/V (q)IIo(w,q) 



(3) 



where Gr/a (p, uj, a) are the retarded/advanced Green 
functions of electrons with spin direction a — ±. Func- 
tion F(u>, a) in the uniform state described by the chem- 
ical potential shifts Sfj, can be written as 



F (uj, a) = 1 — 2n a (u) = tanh 



2T 



(5) 



We see that the integral over u of the first term in ex- 
pression @ proportional to F(e, a) is equal to zero since 
all diffusion poles are in the region Im(u;) < 0. Substi- 
tuting ([5]) into @ we note that expression © is equal 
to the equilibrium one under the substitution of energy 
e — aSfj,/2 by energy e in the expression for correction to 
DOS for spin direction a. Therefore, the total correction 
to DOS for spin direction a has the form 



6v a (e) 



v \ f duj 
— Im > / — 



2- 



C/fl iCt (w,q) tanh ( - 



2T 



[Dq 2 — iuj — ia(tl z + £l p )] 



1 (U Rfi (u, q) - V R (u), q)) tanh (- 

2 [D<f - iui] 2 



-uj — aSfi/2 
2T 



(6) 



where we define 

Ur,m(u>, q 



A 



1 - \vU M (uj,q) 
and M = (0,a), while the polarization operator is 

Dq 2 - iMCLp/Xu 



n M (w,q) 



Dq 2 -icj- iM(fl z + Q p ) 



(7) 



(8) 



Here we include Zeeman splitting il z and exchange en- 
ergy that is self-consistently defined as 



A f de 



(9) 

We assume uniform spin accumulation and magnetic field 
applied parallel to the non-equilibrium magnetization. 

I 



Magnetic field polarizes electron spin in the system and 
as a result the spin polarization S in the paramagnet is a 
sum of two contributions coming from spin accumulation 
and applied magnetic field. 



S = 



-Pz+Sfj] 



(10) 



1 — \v 

Note, that we always assume S/j, to be positive. At the 
same time f2 z can be either positive or negative depend- 
ing on direction of the magnetic field. 

Let us now consider spin accumulation in the param- 
agnet in the two-dimensional limit. We take into ac- 
count that V^ 1 J2q I (2tt) 2 ' Coulomb interaction 
is Vb(q) = e 2 /q and we also assume small parameter 
Xu < 1. We find that the energy dependence of the 
interaction correction to the density of states in two di- 
mensional conductor at |e|,<5/i > T has the form 



<W,4-( e ) 



4tt 2 D V 4 



D 2 K 4 



lnfle T Sn/2\ r) + — ln(|e qp 6n/2\ r) + Xu lnfle ± (0, 



fi p + (5/i/2)|- 



(11) 
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where k = 2ire 2 v is the inverse static screening length. 
Comparing this expression with that in the equilibrium 



case [141 ] . we find that when spin polarization S is zero, 
non-equilibrium spin accumulation splits the singulari- 
ties of the density of states at e = [14] to e = ±Sfi/2 
for electron with spin up and down, correspondingly. If 
S 7^ then the density of states exhibits singularities 
located at e = ±|<J/z/2 -I- Sl z \ for small \v in addition to 
the singularities at ±5/x/2. This case is similar to the 
equilibrium one in the presence of magnetic field. 



CONDUCTIVITY 
Metallic region 

When considering the interaction correction to conduc- 
tivity in the presence of chemical potential shifts we can 



use the results obtained for equilibrium 1J, ll5j similarly 
as it was done for DOS. The combined effect of exter- 
nal magnetic field and spin accumulation on conductiv- 
ity is determined by total spin polarization. Expression 
for conductivity contains function ■Jh[u> coth(a»/2T)] that 
appears as a result of integration of the equilibrium Fermi 
distribution functions. Instead of these equilibrium func- 
tions we obtain for S ^(-dependent part of conductivity 



/oo 
de [F (e, a) F (e + id, -a) - 1] = 
-OO 

(^) (-) 



1 d 

Yduj 



As a result, we obtain the expression for the spin accu- 
mulation dependent interaction correction to the conduc- 
tivity 



8(7 = — i——^-—— Do 2 I did- 

ird V ^ 1 

q,a— ± 



^fl,a(w, q)<9 w (ui — aSfi) coth 



2T 



(Dq 2 — iid — ia(Cl z + £l p )Y 



(13) 



Expression (113[) coincides with that in the equilibrium 
case (l5j with the corresponding spin polarization 5, sim- 
ilarly as it was done for DOS . The interaction correction 
to the conductivity of two-dimensional metallic system 
in the limit of large spin accumulation 8fi > T becomes 



o(5n) - cr(0) = - 



\ve 2 
Air 2 h 



In |#/xt| 



(14) 



while in the small spin accumulation S/i < T regime the 
interaction correction takes the form 



a(Sfi) - <t(0) = -0.08 



Xve 2 ( 5/i 



4ir 2 h V T 



(15) 



If the external magnetic field is applied to the system 
then Spin expressions (TT4|) and (fTB"]) must be substituted 
with |(5/i + il z |. However, when considering the frequency 
dependent conductivity the equivalence between S/x and 
Zeeman energy will be lost. 



Hopping region 

Let us consider hopping conductivity in the presence of 
spin accumulation. Spin dependent contribution to the 
hopping conductivity arises if electron hopping involves 
states that permit double onsite occupancy in the pres- 
ence of Coulomb repulsion 

HE E3|- According to [IB, E3] 
hopping rates are determined by correlation functions of 
occupation numbers of the sites involved in hopping. The 



one site density matrix is given as 
ei(n + +rt_) Un + n- 



exp 



T 



T 



5[i(n + — rt_) 
2T 



(16) 

Here is the one electron energy of localized state i, 
counted from the Fermi level and U > is the intra-site 
electron repulsion potential. 

Expression (|16p coincides with equilibrium one for elec- 
trons in magnetic field. After mapping the approach de- 
veloped in IE 13 to the case of spin accumulation the 
hopping rate from site i to site k can be written as sum 
of contributions proportional to 

1. exp(— [ei ± Sn/2]/T)/ZiZk if site i has one electron 
and site k is empty, 

2. exp(— [ej + ek]/T)/ZiZk if both sites are single oc- 
cupied, 

3. exp(— [2ei + U\/T) / ZiZk if site i has two electrons 
and site k is empty, 

4. exp(-[2e. i + e k + U T 8^/2}/T)/Z i Z k if site i has 
two electrons and site k has one electron. 

where partition function in the case of chemical potential 
shifts is given as 



Zi 



1 



2 cosh I — I < 
2T 



■ei/T 



-(2e 4 +!7)/T 



(17) 



for U > assumed to be equal for all sites. 
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Scaling expression for the variable range hopping con- 
ductivity in this problem was obtained in 18]. For small 
values of spin accumulation defined by the inequality 
< !> where T is the characteristic Mott 
temperature and d is the sample dimension, we obtain 
the expression for the spin accumulation 5(1 dependent 
part of conductivity 
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In 



a(Sfx) 



a{0) 



9s9d 
9S + 9D T 



(18) 



where g$ and gn are the densities of states with energies 
e a and e a + U near the Fermi level correspondingly. Note 
that in this regime In \a(Sfi)/a(Q)\ does not depend on 
dimension d. Again, expression (fT5|) has to be modified 
in the presence of magnetic field. Spin accumulation term 
S/j, must be substituted with |<5/i + fl z \. 



CONCLUSIONS 

The singularities of DOS can be measured with a tun- 
nelling probe electrode, where the singularities can re- 
veal in the zero bias dips of the tunnelling conductivity 
14| . We note that in spin accumulation regime both sides 
of the tunnelling contact will acquire chemical potential 
shifts. In order to resolve the splitting of singularities in 
DOS, the chemical potential splitting in the probe elec- 
trode has to be much smaller than that in the studied 
system. This can by achieved if one takes the probe elec- 
trode with much stronger spin relaxation compared to 
the spin relaxation in the studied system. 

Let us now discuss the experimental observability of 
considered correlation effects. An important issue is the 
value of chemical potential shift 5/j, compared to the value 
of the inverse spin relaxation time 1/t s in the system. 
The applicability of expressions (fTTj). (fT4|) and ([T5|) re- 
quires \J5ht s to be large. We take the spin relaxation 
time to be t s ~ 100 ps and the spin accumulation Sfi ~ 1 
meV. The 1 meV spin accumulation requires correspond- 
ingly low < 10 K temperatures. We obtain \f5JvF s ~ 10. 
Also assuming the diffusion coefficient to be 10 cm 2 s _1 , 
we estimate l s = \/Dt s ~ 300 nm and yjD/\8^\ ~ 30 
nm. The above conditions on the magnitude of the spin 
relaxation length and spin accumulation are realizable ex- 
perimentally in the ferromagnet - doped semiconductors 
contacts 



the nature of this l arg e 
value of spin accumulation is still under debate 22J, |23j • 
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